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Abstract 

The first-order cosmic ray acceleration at ultrarelativistic shocks is investigated us¬ 
ing the Monte Carlo method. We apply a method of discrete particle momentum 
scattering as a model of particle pitch angle diffusion to reproduce highly anisotropic 
conditions at the shock wave. Shocks with Lorentz factors 7 up to 320 and varying 
magnetic held inclinations ip are considered. Values of diffusion coefficients upstream 
in the point where energy spectral indices stabilize to the limit 2.2 were calculated. 
The obtained acceleration time does not depend on shock conditions. 

Introduction. 

Till quite recently the discussion of hrst order Fermi acceleration at relativistic 
shocks was restricted to small Lorentz 7-factors (for review see Ostrowski[l] and 
Kirk[2]). Bednarz and Ostrowski[3] found a simple mechanism that allows particles 
to be accelerated at shock waves with large Lorentz factors. It implies that an incli¬ 
nation of formed spectrum does not depend on downstream but only on upstream 
conditions and a shock velocity. The mechanism is a natural proccess accelerating 
particles in astronomical objects where such shocks occur. 

Among these phenomena where ultrarelativistic shocks are anticipated are gamma- 
ray bursts. They are non-thermal bursts of low energy 7-rays. Typical GRB lasts 
for about 10 sec. Their broad spectra usually peak between a few 100 keV and a few 
MeV (for review see Piran [4]). The Lorentz factor of an expanding hreball larger 
then 300 is required in order to the pair production optical depth to be sufficiently 
small (Baring [5]). Recently the host galaxy with redshift 3.42 was indentified for 
GRB971214 (Kulkarni et ah [6]) which confirms the suggestion that they originate 
at cosmological sources (Meegan et ah [7]). 

Numerical simulations of the acceleration process. 

In simulations we follow the procedure used in Bednarz and Ostrowski [8]. Seed 
particles are injected downstream the shock with the same initial energy. Theirs 
trajectories are derived in homogenous magnetic held Bq perturbed by inhomogeni- 
ties which are simulated by particle momentum scattering within a cone with small 
angular opening AG less than the particle anisotropy ~ I/7 (see Ostrowski [9]). A 
particle is excluded from the simulation either if it escapes through the free-escape 
boundary placed far downstream the shock or if it reaches a time or energy larger 
than the assumed upper limit. These particles are replaced with ones arising from 
splitting the remaining high-weight particles with preserving their physical param¬ 
eters. They are not replaced if they reach upper limit of time. 
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Figure 1: The simulated spectral indices a for particles accelerated at shocks with 
different upstream diffusion represented by logio(ft;±/ft;||)- Downstream we assume 
K,± = 0. Perpendicular dashed line points to ‘boundary point’. 



Figure 2: The simulated acceleration times tacc for particles accelerated at shocks 
with different Lorentz factors 7 . The considered upstream magnetic held inclination 
are represent by different lines with indicated simulation errors. The tacc value near 
1.5 is for 'ijj = 30° and results from large simulation error at 7 = 20. 





Computations are performed in the respective - upstream or downstream - plasma 
rest frame. If particles cross the shock their parameters are transformed to the shock 
rest frame. The respective contribution divided by their velocity component normal 
to the shock added to momentum bin forms particle spectra with shifting cut-off 
point toward higher energies. These spectra allow us to derive their inclination and 
acceleration time (see Bednarz and Ostrowski [ 8 ]). Time measured in the shock rest 
frame in downstream r^/c units was transformed to the downstream plasma rest 
frame {vg is a particle gyroradius). That can be done because we know the ratio of 
mean time that particle spends upstream to mean time it spends downstream. This 
value is small and due to this the time measured in the shock rest frame is of 1.05 
factor shorter than downstream for large 7 . 

In this paper we simulate different conditions upstream and downstream the 
shock, represented by cross-held diffusion coefficient k± and the parallel diffusion 
coefficient k,\\. 

Results. 

We take into consideration shocks with the Lorentz factors equal 7 =20, 40, 80, 
160, 320 and magnetic held inclinations upstream the shock ip = 15°, 30°, 45°, 
60°, 75°, 90°. 

At hrst we consider downstream conditions without magnetic huctuations. We 
search by simple data inspection as in Fig. 1 for a value where the spectral 

index reaches limiting value of 2.2. It yields a ht: 

^±/^|| = 0.25 ■ 7 “^-^^'^ 

The ht is good in full considered parameter space. 

The acceleration time tacc for the boundary points is presented in Fig. 2. where 
one can seen the lack of any change with ip. The hgure exibits also a slow tacc decrease 
from 1.2 to approach the value 1.1 at 7 = 80. We observe small tendency of tacc 
to grow if a increases to 2.3-2.4 and no further change if kj_/k\\ grows. In the next 
step we introduce magnetic held huctuations downstream to be k,_l/k,\\ = 7.1 ■ 10“"^ 
or 1.1 ■ 10 “^. tacc behaves nearly as in the case without huctuations, with only a 
small reduction (circa 0.1) for the larger huctuation. Thus the resulting acceleration 
times are always approximately tacc = 1 (in units of r^/c downstream the shock) 

Discussion. 

We have obtained the result that is not suprising if we recall the way particles are 
accelerated (Bednarz and Ostrowski [3]) and that they do not spend much time 
upstream. 

We expect our results can be applied in modelling relativistic shock waves ex¬ 
pected to occur in gamma-ray burst sources and at pulsar wind terminal shocks. 
The model of Kenel and Coroniti [10] assumes for Crab Nebula initial how with 
Lorentz factor equal to 10®. One should note that the acceleration time is short, 
close to the value of r^/c downstream the shock. However, it is not very short with 
energy gains ~ 7 ^ at single particle-shock interaction. It results from the fact that 
it is not possible to rehect a particle from the ultrarelativistic shock ( cf., Begelman 
and Kirk [11], Ostrowski [1]) 
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